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Abstract. We present a detailed analysis of chemical abundances in a sub-damped Ly a absorber (sub-DLA) at 
z = 2.187 towards the quasar HE 0001-2340 (z cm = 2.28). Our study is based on high-resolution (R m 45,000) 
spectroscopic data from the UVES instrument installed on the ESO Very Large Telescope (VLT). This sub-DLA 
system consists of at least 25 individual subcomponents spanning a restframe velocity range of ~ 400 kms . 
The total neutral hydrogen column density is log AT (Hi) ~ 19.7. Detected species include Cn, Civ, Nl, Nil, Ol, 
Mgn, Alii, Aim, Sin, Siiv, Pn, Fen, and possibly Dl. For the dominating neutral gas component at v TC \ = +13 
kms -1 (relative to z = 2.187) we derive an oxygen abundance of [0/H]= —1.81 ± 0.07 (1/65 solar). With its 
, extremely low nitrogen content ([N/H] < —3.3 and [N/0]< —1.5) the absorber exhibits a classic massive star 

C^*) • abundance pattern. Our measurements place the z = 2.187 absorber towards HE 0001— 2340 among the systems 

with the lowest ever measured [N/a] ratios in the Universe. The low [N/O] value is consistent with the idea that 
primary nitrogen production by the very first stars have enriched the intergalactic gas to a level of [N/0]~ —1.5. 
Peculiar abundances are found in the outermost blue components near —350 kms 1 (in the z = 2.187 restframe) 
where we observe significant overabundances of phosphorus ([P/C]~ +1-5), silicon ([Si/C]~ +0.6), and aluminum 
([Al/C]~ +0.5) after correcting for the effects of ionization. Our study suggests that the sight line passes through 
|' the gaseous environment of one or more stellar clusters that locally enriched their interstellar neighbourhood by 

1~ \ supernova ejecta generating the observed abundance anomalies. The large velocity spread of the entire absorption 

system points to a merger that triggers the formation of these clusters. We thus may be observing a young galaxy 
at z m 2 that currently forms out of a merger event. 
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5h I 1. Introduction the chemical enrichment by the first generations of stars 

■ - '< (e.g., Lu et al. 1996; Prochaska & Wolfe 1999). Such ab- 

Absorption line systems at high redshift provide a wealth SQrbers therefore may represent progenitors of present-day 

of information about the chemical evolution of the early g a l ax ies 
Universe. Next to the Ly a forest and the Lyman limit sys- 

terns that trace the low-density IGM or the intergalactic Most of the information about abundances of metals, 

environment of galaxies at neutral hydrogen column den- dust ' and molecules m these systems has been obtained 

sities of JV(Hi) < 10 19 cm" 2 , the physical and chemical from the dam P ed ^ a systems (DLAs) , which have Hi 

properties of absorbers with N(B i) > 10 19 cm ~ 2 are par- column densities > 2 x 10 cm (e.g., Ledoux, Petitjean, 

ticularly interesting. These systems contain most of the & Snanand 2003). Absorption systems with somewhat 

neutral gas mass in the Universe at z > 1 (Lanzetta, lower column densities (10 19 cm < AT(Hl) < 2 x 10 20 

Wolfe, & Turnshek 1995; Wolfe et al. 1995) and consist cm ) are the so-called sub-damped Lya systems (sub- 

of heavy elements in an abundance pattern that suggests DLAs ' f*' Dessauges-Zavadsky et al.2003; Peroux et 

al. 2003). While there is no compelling evidence that DLAs 

~~Z ', ~Z ~~ T ~ r, t-, . , , and sub-DLAs trace different classes of protogalactic ob- 

bend offprint requests to: P. Hicnter r ° 

e-mail: prichter@astro.uni-bonn.de J ects ' abundance measurements in sub-DLAs are of par- 

* Based on observations carried out at the European ticular interest because of their absorption characteristics. 

Southern Observatory (ESO), La Silla, under prog. ID No. For total H I column densities < 2 x 10 20 cm ~ 2 and metal- 

166.A-0106 with the UVES spectrograph at the ESO Very licities less than ~ 0.1 solar, the absorption lines from neu- 

Large Telescope, Paranal, Chile. tral oxygen and neutral nitrogen (two important elements 
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Fig. 1. Portion of the UVES spectrum of HE 0001-2340 in the wavelength range between 3760 and 3912 A. Absorption 
by Sin AA1190.4, 1193.3, SimA1206.5 (blended by intervening H i), and Hi Lya A1215.7 from the z = 2.187 sub-DLA 
is shown in the upper, middle, and lower panel, respectively. 



for disentangling the complex chemical enrichment history 
of high-z absorbers) are relatively weak. This enables us 
to derive accurate abundances for these two important el- 
ements in sub-DLAs without facing the problem of heavy 
line saturation in the oxygen lines. 



In this paper we present the detailed analysis of a sub- 
damped Lya system at z = 2.187 in the direction of the 
quasar HE 0001— 2340. This paper represents the first one 
in a series of two dealing with the chemical and physi- 
cal properties of particularly interesting sub-DLA systems 
in the data set of the ESO-VLT Large Programme "The 
Cosmic Evolution of the IGM". This paper is organized 
as follows: in Sect. 2 we present the observations, the data 
reduction, and the data analysis method. Absorption by 
neutral and weakly ionized species is presented in Sect. 3. 
Sect. 4 deals with the properties of the highly ionized gas. 
We discuss the results in Sect. 5 and summarize our study 
in Sect. 6. 



2. Observations, data handling, and analysis 
method 

The quasar HE 0001-2340 (V = 16.7, z cm = 2.28) was 
observed with the UVES instrument installed on the VLT 
as part of the ESO-VLT Large Programme "The Cosmic 
Evolution of the IGM" , which aims at providing a homo- 
geneous data set of high- resolution (R ~ 45,000), high 
signal-to-noise (S/N> 30) spectral data for QSO absorp- 
tion line spectroscopy of the IGM in the rcdshift range 
between z = 2.2 — 4.5. 

HE 0001-2340 was observed through a 1" slit (with a 
seeing typically < 0.8") with two setups using dichroic 
beams splitters for the blue and the red arm (Dicl, 
B346nm+R580nm, and Die 2, B437nm+860nm, respec- 
tively). This setup provides a wavelength coverage from 
- 3050 to 10,400 A, with small gaps near 5750 and 
8550 A. The total integration time (Die 1 + Die 2) for 
HE 0001— 2340 was 12 hours. The raw data were reduced 
using the UVES pipeline implemented in the ESO-MIDAS 
software package. The pipeline reduction includes flat- 
fielding, bias- and sky-subtraction, and a relative wave- 
length calibration. The individual spectra then have been 
coadded and corrected to vacuum wavelengths. The S/N 
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Fig. 2. Absorption profiles of Cn and Alii from the 
z = 2.187 sub-DLA are displayed (velocities refer to the 
z = 2.187 restframe). The Cn A1334.5 absorption (up- 
per panel) demonstrates the complex velocity structure 
of this absorber. 16 individual components are identified, 
given by the letters A to P. These components are indi- 
cated by the vertical dashed lines. The lower panel shows 
an example for a multi-component Voigt profile fit from 
the FITLYMAN routine. Absorption by Alii A1670.8 (dots) 
is overlayed by the corresponding optimum fit (solid line). 
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3. Neutral and weakly ionized species 

3.1. Velocity structure and b values 

Fig. 2 (upper panel) shows the Cn A1334.5 line plotted 
on the z = 2.187 restframe velocity scale, demonstrating 
the complex component structure of this absorption sys- 
tem. The Cn line, among others, is sufficiently strong to 
show also weak velocity components and lies in a region 
that is not affected by line blending. Moreover, this line 
has high S/N, so that the various subcomponents can be 
reliably identified. To derive a unique model for the ve- 
locity component structure we have fitted Voigt profiles 
to a number of high S/N lines of Cn, Mgn, and Sin. We 
identify 16 velocity subcomponents, in Fig. 2 shown with 
vertical dashed lines. In the z — 2.187 restframe these ab- 
sorption components have relative velocities of +13, —29, 
-111, -145, -159, -187, -217, -229, -240, -271, -284, 
-305, -327, -340, -353, and -369 kms" 1 . This velocity 
model then has been used to derive b values and column 
densities for each component via Voigt-profile fitting. The 
resulting values for b and log TV and their la errors are 
listed in Table 1. Throughout the following, we label the 
individual velocity components with the letters A to P in 
alphabetical order (i.e, component A is the one at +13 
kms -1 ). The strongest absorption in Cll and other neu- 
tral and weakly ionized species is seen in component A, 
thus at the outermost red wing of the overall absorption 
structure. 



in the spectrum generally is very high and attains a max- 
imum of ~ 150 per resolution element for A > 4000 A. In 
Fig. 1, a representative portion of the UVES spectrum of 
HE 0001—2340 in the wavelength range between 3760 and 
3912 A is displayed, showing selected absorption lines of 
the sub-DLA system at z = 2.187. 

For the spectral analysis of the data we made use of 
the FITLYMAN package in MIDAS (Fontana & Ballester 
1995). This routine uses a \ 2 minimization algorithm to 
derive column densities, N, and Doppler-parameters, 6, 
via Voigt-profilc fitting, taking into account the spectral 
resolution of the instrument (see above). First of all, we 
have disentangled the component structure in the z = 
2.187 absorber for a) neutral and weakly ionized species 
(singly ionized species plus Aim), and b) highly ionized 
species (Civ and Siiv), as explained in detail in Sects. 3 
and 4. Using a number of high S/N absorption profiles 
from various species we then obtained Doppler parameters 
for each component by Voigt profile fitting. We use a single 
b value per absorption component to fit all neutral and 
weakly ionized species except hydrogen and deuterium. 
Fig. 2, lower panel, shows an example for a Voigt profile fit 
to the Al II A1670.8 line. The spectral data are plotted with 
black points, the corresponding fit is shown as solid line. 
Throughout this paper we use wavelengths and oscillator 
strengths from the atomic line list of Morton (2003). 



3.2. Metals 

A number of neutral and weakly ionized species are de- 
tected in the various velocity components, including Cn, 
Ni, Nil, Oi, Mgn, Al ii, Aim, Sin, Pn, and Fen. A se- 
lection of velocity profiles of these species are shown in 
Fig. 3. Column densities for these ions as well as upper 
limits for undetected species (e.g., Znii, Sil, and molecu- 
lar hydrogen, H2) are listed in Table 1. The strong lines 
of Cn (AA1334.5, 1036.3) and Mgn (AA2803. 5, 2796.4) are 
saturated in the dominating neutral components A, E, K, 
O, and P. As there are no weaker transitions available for 
these species, the column densities derived for C 11 and 
Mg 11 in these components are relatively uncertain (see 
Table 1). For Oi, Ni, Alii, Aim, Sin, and Fen unsatu- 
rated lines are available that allow us to derive accurate 
column densities from profile fitting. N 1 absorption is de- 
tected solely in component A. The measured equivalent 
width of the Ni feature in the A1199.6 line (see Fig. 3, 
third panel in the left row) is only 1.8 mA, correspond- 
ing to a 3.0cr significance at the measured local S/N of 
<~ 130 per resolution element. Remarkable is the detec- 
tion of Pn in outermost blue components O and P, as 
shown in Fig. 4. This is surprising, since P 11 has a rela- 
tively low cosmic abundance (log (P/H)q = —6.44) and 
is not observed in the strongest neutral component A. 
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Fig. 3. Absorption profiles of Cn, Ni, Nil, Oi, Mgn, Alii, Aim, Sin, Pn, S II, Fen, and Fein for the z = 2.187 
sub-DLA are shown. The individual absorption components are indicated with the vertical dashed lines. The velocity 
scale refers to the z = 2.187 restframe. 
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Fig. 4. Zoomed-in version of P n A1152.8 absorption in 
the z — 2.187 absorber in the range between —450 and 
—250 kms" 1 . Pn absorption in components O and P (in- 
dicated by the dashed lines) is clearly detected. The solid 
line shows the optimum two-component Voigt-profile fit. 



3.3. Neutral hydrogen absorption 

Neutral hydrogen absorption in the z = 2.187 sub-DLA 
is seen in Lya, Ly/3, and Ly7, as shown in Fig. 5. The 
other lines from the H I Lyman series either have very low 
S/N or are located bluewards of the available UVES wave- 
length range. As clearly visible in Figs. 1 and 5, the Lya 
absorption is very asymmetric due to the large number of 
neutral gas components at the blue side that contribute 
to the overall H I absorption. The red part of the Lya ab- 
sorption shows a well-defined damping wing caused by the 
strongest neutral gas component at +13 kms -1 (compo- 
nent A). Due to the asymmetry, a free single-component 
fit of the Lya, Ly/3, and Ly7 absorption cannot even 
roughly reproduce the observed shape of the absorption. 
Therefore, the complex velocity structure of the z = 2.187 
absorber has to be taken into account for the profile fit 
of the H I Lyman series absorption. Fixing the centroid 
of the fit at +13 kms -1 and only fitting the damping 
wing on the red side of the Lya, Ly/3, and Ly7 absorp- 
tion simultaneously results in an H I column density of log 
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Table 1. Logarithmic ion column densities and b values for the z = 2.187 absorber towards HE 0001-2340 
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12.94 ± 0.10 


3.5 ± 1.3 






Alii 
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Mgii 
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12.15 ± 0.04 
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11.00 ± 0.05 








Sin 
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la errors and 3<r upper limits are given. 



7V(Hi)a = 19.60 ±0.4 for component A. Similarly, we can 
constrain the total Hi column density in the two bluest 
velocity components O and P to log iV(Hi)op < 18 by 
fixing the velocity centroid at —361 kms -1 and fitting the 
blue wing of the Ly j3 and Ly 7 absorption. Note that the 
additional absorption occurring in the blue wing of the 
Lya, Ly/?, and Ly7 absorption is possibly due to deu- 
terium (D 1) in components O and P in combination with 
an additional weak H 1 component near —440 km s _1 . This 
will be further discussed in the Appendix. 

The entire Hi absorption in Lya, Ly/3, and Ly7 pro- 
file can be modeled by taking into account the absorp- 
tion from the intermediate velocity components. For our 
model, we consider absorption from the strongest neutral 
gas components A, B, E, K, L, O, and P, for which we 



have fixed the velocity-component structure as seen in the 
metal-line absorption. Characteristic is the presence of a 
residual flux spike in the Ly/3 profile near —300 kms -1 , 
which needs to be reproduced by our model. A similar 
spike is not seen in the Ly7 absorption profile, but it is 
likely that in this region the Ly7 line is blended with 
absorption from the Lya forest. The best fitting model 
has logarithmic H 1 column densities of 19.60, 17.79, 18.25, 
15.03, 17.71, 17.53, and 17.32 for the components A, B, 
E, K, L, O, & P, respectively. This model is shown Fig. 5 
with a thick solid line; the expected flux spike near —300 
kms -1 in the Ly7 absorption is indicated with a thick 
dashed line. Note that this fit for the intermediate veloc- 
ity components is not unique, as other combinations of 
H 1 column densities also can reproduce the shape of the 
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Fig. 5. Hi Lya, Ly (3, and Lyj absorption in the z = 2.187 absorber. The thick solid line represents the best-fitting 
H i model, as described in Sect. 3.3. 



Table 2. Summary ol logarithmic metal abundances in 
component A 



Species 


LP. 

[eV] 


log (X/H) a 
(+12) 


[X/H] b 


Cn 


24.38 


8.52 


-1.76 ±0.21 


Ni 


14.53 


7.95 


-3.50 ±0.05 


Nil 


29.60 


7.95 


< -2.75 


Oi 


13.62 


8.69 


-1.81 ±0.07 


Mgn 


15.04 


7.58 


-2.05 ±0.11 


Alii 


18.83 


6.49 


-2.11 ±0.06 


Aim 


28.45 


6.49 


-2.64 ±0.06 


Sin 


16.35 


7.56 


-1.81 ±0.07 


Pn 


19.77 


5.56 


< -1.18 


Sn 


23.34 


7.20 


< -1.68 


Fen 


16.19 


7.50 


-2.12 ±0.05 


Fein 


30.65 


7.50 


< -2.23 


Znii 


17.96 


4.67 


< -1.23 



a Reference abundances are taken from Allende Prieto et 
al. (2001) for oxygen and from the compilation of Morton 
(2003) and references therein. 
b [X/H]=log (X/H)-log (X/H) . 



observed Hi absorption. An important outcome of this 
model is, however, that the absorption of the intermedi- 
ate components does not significantly influence the shape 
of the H I absorption in blue and red wing of the Lyman 
lines. Therefore, the values of log AT(Hl) A = 19.60 ± 0.4 
and log iV(Hl)op < 18 derived above represent reliable 
estimates for the Hi column densities in components A 
and O+P. 



3.4. Abundances and abundance ratios 
3.4.1. Component A 

First, we concentrate on the dominating neutral compo- 
nent (component A) in our discussion of abundances and 
abundance ratios in the z = 2.187 absorber. The best ele- 
ment to derive the overall abundance in the gas is oxygen, 
since neutral oxygen and neutral hydrogen have very sim- 
ilar ionization potentials and both elements are coupled 
by a strong charge-exchange reaction. Moreover, oxygen 
does not deplete significantly into dust grains. Our anal- 
ysis yields [O/H] = —1.81 ± 0.07 1 , which corresponds 
to an abundance of ~ 1/65 solar. In striking contrast, 
the nitrogen abundance is significantly lower, [N/H] = 
-3.50 ± 0.05 or - 1/3200 solar, as derived from Ni. 
The resulting nitrogen-to-oxygen ratio [N/O] = —1.69, 
if accurate, would be one of the lowest ever measured 
in DLAs so far. However, the interpretation of Ni/Hl is 
not as straightforward as for Ol/Hl. Although the ioniza- 
tion potential of Ni is only slightly higher than that of 
hydrogen (14.53 versus 13.60 eV), ionization effects can 
significantly affect the determination of [N/H] from Ni 
and H I in absorbers with overall column densities < 10 20 
cm~ 2 . The reason for this is that charge exchange reac- 
tions are less efficient between N I and H I (in contrast 
to Ol and Hi; see Butler & Dalgarno 1979) and Ni has 
a relatively large photoionization cross section (Sofia & 
Jenkins 1998). Thus, a considerable fraction of the nitro- 
gen in component A may reside in the N + phase, partic- 
ularly if a substantial ionizing radiation field is present. 



1 [X/H]=log (X/H)-log (X/H) 
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While the uncertain column density constraint from the 
blended Nil A1084.0 line (log JV(Nll) < 12.80) does not 
provide a very stringent estimate of the abundance of 
singly ionized nitrogen, we can check upon possible ion- 
ization effects in component A using other weakly ionized 
species. The ratio of (Aim/An) w 0.3 suggests that ion- 
ization could be important. This ratio, however, is difficult 
to interpret due to the uncertain recombination rate of 
A1+ (sec Vladilo et al. 2001; Nussbaumer & Storey 1986) 
and the observed discrepancies between (Aim/ An) and 
other ionic ratios in high-z photo-ionization models (e.g., 
D'Odorico & Petitjean 2001). An important element to 
investigate ionization conditions in interstellar and inter- 
galactic gas is the a element sulfur (Richter et al. 2001). 
Sn has an ionization potential of 23.23 eV and depletes 
only very little into dust, so that we expect most S to 
be in form of Sn. Since the ionization potential of Sn 
is significantly higher than that of Oi and Hi (~ 13.6 
eV for both), Sn lives in gas where part of the hydrogen 
is ionized. From the observed Sn/Oi ratio we thus can 
learn about the hydrogen ionization fraction (see Richter 
et al. 2001 for details). Although S n is not detected in our 
sub-DLA, the 3a upper limit for the S n column density in 
component A (log iV(Sll) < 13.12) yields [S/H]< -1.68 
and [Sii/0i]< +0.13. If one assumes a solar sulfur-to- 
oxygen ratio, the possible small excess of S II compared to 
O I implies that the degree of ionization in component A 
is low (< 30 percent). 2 We conclude that Nii/Ni most 
likely is small, so that we can give a conservative limit for 
the nitrogen abundance in component A of [N/H]< —3.3 
and [N/0]< -1.5. 

The apparently low degree of ionization in combination 
with the equal abundances of Si and O speaks against any 
significant dust depletion of silicon in component A. The 
situation is somewhat less clear for the elements Mg, Al, 
and Fe, which are all underabundant compared to oxygen 
([Mg/0]= -0.24, [Al/0]= -0.30, and [Fe/0]= -0.31). 
These elements could be depleted into dust grains. In fact, 
with [Fe/0]= -0.31 and [Si/0]= 0.00 the absorber follows 
the trend for dust depletion patterns found in damped 
Lya systems and in the Galactic halo (see Petitjean, 
Srianand, & Ledoux 2002, their Fig. 7). However, next to 
dust depletion nucleosynthesis effects certainly have to be 
considered for these elements given the low overall abun- 
dance and the very low nitrogen content. In this context, 
the observed abundances of iron and aluminum could also 
be explained by a slightly enhanced [a/Fe] element ra- 
tio and a mild odd-even effect for Al. In contrast to Fe 
and Al, the column density of magnesium unfortunately 
is too uncertain (due to saturation of the two available 
Mg II lines) to draw meaningful conclusions about the ap- 
parent Mg underabundance. The absence of molecular hy- 
drogen absorption in component A represents yet another 



2 Note that if the intrinsic S/O ratio in this system is sig- 
nificantly lower than solar, the ionization fraction would be 
underestimated by this method. Some DLAs indeed show sub- 
solar S/O ratios (e.g., Prochaska, Howk, & Wolfe (2003)). 



piece of information on the dust abundance in the gas. 
From the H 2 column density limit (see Table 1) we ob- 
tain an upper limit for the molecular hydrogen fraction 
of / = [2JV(H 2 )/(iV(Hi) + 2AT(H 2 ))] » 5 x lO" 7 . As H 2 
preferentially forms on the surfaces of dust grains, the 
lack of H 2 further supports the above conclusion that the 
dust abundance in component A must be small (see also 
Ledoux, Petitjean, & Srianand 2003). Little can be said 
about the carbon abundance in component A, as the value 
of [C/H] = —1.76 ± 0.21 - although similar to that of oxy- 
gen - is very uncertain due to saturation of the two C n 
lines. If the b value would be only slightly higher than the 
assumed 4.5 kms -1 , the Cn column density (and thus the 
derived carbon abundance) would drop substantially. 

All measured abundances in component A are summa- 
rized in Table 2. 

3.4.2. Components O and P 

We now turn to components O and P, which show a very 
distinct abundance pattern. As discussed in Sect. 3.3, the 
Lya Voigt profile fit yields log iV(Hi) < 18 for compo- 
nents O and P together. The combined Oi column den- 
sity for O and P is log N « 12.7 (see Table 1), from 
which we infer an oxygen abundance of [0/H]op] > —2.0. 
This is consistent with the oxygen abundance derived for 
component A. The large value for log (Cii/0i)= +1.15 
in component O implies that the degree of ionization 
must be large. Note that the Cn absorption in compo- 
nents O and P is not saturated (unlike in component 
A) , so that the derived C n column densities are accu- 
rate (see Table 1). From Cn, Mgn, Sin, Alii, and Pn 
we obtain (Mgii/Cn)op = -0.76, (Siii/Cn) p = -0.54, 
(A1ii/Cii)op = -1.21, and (Pii/Cii) p = -1.45. From 
Fen and Fern we derive ((Feii+Feiii)/Cn)op = -0.74. 
These ratios suggest significant overabundances of sili- 
con, aluminum, and phosphorus in components O and 
P. However, a confirmation of this statement and a more 
quantitative estimate of the apparent overabundances has 
to await a detailed photoionization modeling for these 
species. This will be presented in Sect. 4.2. 



4. Highly-ionized species 

4.1. Component structure and column densities 

Absorption by highly ionized gas in the z = 2.187 absorber 
is seen in C iv AA1550.8, 1548.2 and Si iv AA1402.8, 1393.3. 
Ovi and Nv are not detected, but blending problems 
hamper a reliable detection of these ions. The compo- 
nent structure seen in C iv and Si IV (Fig. 6) is even more 
complex than that of the neutral and weakly ionized gas. 
All in all, we identify 25 absorption components including 
the 16 components detected in neutral and weakly ionized 
species. All C iv and Si iv measurements are summarized 
in Table 3. 
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Fig. 6. Absorption by Civ (AA1550. 8, 1548.2) and Siiv 
(AA1402.8, 1393.8) in the z = 2.187 absorber is shown. 
The absorption pattern of these high ions is even more 
complex than for the neutral and weakly ionized species. 
We identify 25 components that partly correspond to the 
absorption components seen in Cn (see also Table 3 for 
details). 



4.2. Photoionization modeling for comp. O and P 

Our measurements of C n, Civ, Sin, and Siiv allow us 
to model the photoionization in components O and P and 
to constrain the ionization parameter U = n 7 /nn, the 
ratio of ionizing photon density to gas density. We have 
used the photoionization code CLOUDY v94.00 (Fcrland 
1996) and have modeled a plane-parallel slab of gas with a 
mctallicity of 0.01 solar and solar relative chemical abun- 
dances (Fig. 7). The calculations were done using a con- 
stant density and were stopped whenever the neutral hy- 
drogen column density reached 10 18 cm~ 2 . The total hy- 
drogen column density in the model, 7V(H i+H n), comes 
out to 1.0 x 10 20 cm~ 2 . For the ionizing radiation, we used 
the Madau, Haardt & Rees (1999) metagalactic spectrum 
at z = 2.2. Our CLOUDY modeling for the column den- 
sity ratios of C iv/C n and Si iv/Si n implies ionization pa- 
rameters of log U = —2.59 and —2.35, respectively, thus 
in good agreement with each other (see Fig. 7, vertical 
dashed lines). A value of log U « —2.5 further is consis- 
tent with the observed column densities/column density 
limits for Fell and Fein. From the Alin/Aln ratio we ob- 



Comp. 


"rel 

[kms -1 ] 




Ion 


log AT" 


6 a 

[kms -1 ] 


A 


+ 13 


2.18714 


Civ 
Siiv 


12.01 ± 0.05 
12.28 ± 0.04 


4.5 ± 0.7 

5.6 ± 0.4 


Al 


+5 


2.18705 


Si iv 


13.01 ± 0.03 


9.4 ± 1.1 


A2 





2.18700 


Civ 


13.58 ± 0.04 


9.7 ± 1.4 


A3 


-11 


2.18688 


Siiv 


12.52 ± 0.08 


18.5 ± 3.6 


B 


-29 


2.19669 


Civ 


13.48 ± 0.06 


14.6 ± 2.7 


Bl 


-53 


2.18644 


Civ 
Siiv 


12.84 ± 0.11 
12.11 ± 0.14 


11.6 ± 4.4 
3.5 ± 1.8 


B2 


-76 


2.18619 


Civ 


13.02 ± 0.12 


8.8 ± 4.9 


B3 


-88 


2.18606 


Civ 
Siiv 


13.22 ± 0.07 
12.45 ± 0.04 


6.8 ± 0.9 
4.7 ± 0.3 


C 


-111 


2.18582 


Civ 
Siiv 


12.93 ± 0.06 
11.85 ± 0.05 


6.5 ± 0.8 
4.9 ± 1.3 


CI 


-131 


2.18561 


Civ 


12.63 ± 0.06 


7.6 ± 1.1 


D 


-145 


2.18546 


Civ 
Siiv 


12.40 ± 0.11 
11.74 ± 0.06 


6.5 ± 2.4 
6.4 ± 1.8 


E 


-159 


2.18531 


Siiv 


11.86 ± 0.23 


15.4 ± 4.4 


El 


-174 


2.18515 


Civ 


11.65 ± 0.07 


9.7 ± 2.0 


F 


-187 


2.18501 


Civ 
Siiv 


12.72 ± 0.10 
11.90 ± 0.06 


11.2 ± 2.1 
11.4 ± 1.7 


G 


-217 


2.18469 


Civ 
Siiv 


12.75 ± 0.04 
12.06 ± 0.05 


8.8 ± 1.8 
3.1 ± 1.1 


H 


-229 


2.18457 


Civ 
Siiv 


12.90 ± 0.04 
12.49 ± 0.03 


9.0 ± 1.3 
6.6 ± 0.5 


I 


-240 


2.18445 


Civ 
Siiv 


12.58 ± 0.08 
11.80 ± 0.09 


9.4 ± 2.2 
9.3 ± 1.5 


J 


-271 


2.18412 


Siiv 


11.68 ± 0.13 


9.8 ± 3.4 


K 


-284 


2.18398 


Siiv 


11.55 ± 0.17 


2.2 ± 1.6 


L 


-305 


2.18376 


Siiv 


11.80 ± 0.10 


5.0 ± 2.8 


LI 


-324 


2.18356 


Siiv 


12.11 ± 0.08 


4.1 ± 1.7 


M 


-327 


2.18352 


Civ 
Siiv 


12.63 ± 0.05 
12.40 ± 0.07 


14.2 ± 3.0 
9.0 ± 1.3 


N 


-340 


2.18339 


Civ 
Siiv 


12.74 ± 0.05 
12.09 ± 0.04 


9.5 ± 2.3 
2.0 ± 0.3 


O 


-353 


2.18325 


Civ 
Siiv 


12.41 ± 0.07 
12.69 ± 0.05 


8.9 ± 2.4 
7.2 ± 2.2 


P 


-369 


2.18308 


Civ 
Siiv 


12.88 ± 0.06 
13.20 ± 0.04 


8.5 ± 1.6 
5.5 ± 1.0 



a 1(T errors and 3cr upper limits arc given. 

tain a significantly higher value of log U = —1.92 (Fig. 7, 
vertical dotted line) . This high value most likely is a result 
of the poorly known dielectric recombination rate of Al + 
at high temperatures (Petitjean, Rauch, & Carswell 1994; 
Nussbaumer & Storey 1986) and therefore will not be con- 
sidered here any further. Comparing the column densities 
predicted by the CLOUDY model with those measured al- 
lows us to constrain the metal abundances in components 
O and P. The results are summarized in Table 4. For oxy- 
gen, carbon, magnesium, and iron we obtain abundances 
that are consistent with those measured in component A. 
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Fig. 7. CLOUDY photoionization model for several ions 
detected in components O and P. The model assumes an 
overall abundance of 0.01 solar together with a Madau, 
Haardt, & Rees (1999) metagalactic ionizing spectrum for 
z = 2.2. The dashed vertical lines indicate the fits for the 
Civ/Cll (left) and Siiv/Sill (right) column density ra- 
tios, the dotted vertical line indicates the fit for Alni/Aln 
(see Sect. 4.2 for details). 

Table 4. Summary of logarithmic metal abundances in 
components O & P 



Species 


log (X/H) a 
(+12) 


[x/nr 


[X/C] b 


C 


8.52 


-1.9 ±0.1 




O 


8.69 


-2.0 ±0.1 


-0.1 ±0.1 


Mg 


7.58 


-1.7 ± 0.1 


+0.2 ±0.1 


Al 


6.49 


-1.4 ±0.1 


+0.5 + 0.2 


Si 


7.56 


-1.3 ±0.1 


+0.6 + 0.2 


P 


5.56 


-0.4 ±0.1 


+ 1.5 + 0.1 


Fe 


7.50 


< -1.9 


< 0.0 



a Reference abundances are taken from Allende Prieto et 
al. (2001) for oxygen and from the compilation of Morton 
(2003) and references therein. 

b Based on the CLOUDY photoionization model described in 
Sect. 4.2. 

For aluminum, silicon, and in particular for phosphorus, 
the CLOUDY model results in relative abundances that 
are significantly enhanced compared to solar. The over- 
abundances relative to carbon amount to 0.5 ±0.1 dex for 
aluminum, 0.6 ± 0.1 dex for silicon, and 1.5 ± 0.1 dex for 
phosphorus (see Table 4). 

5. Discussion 

Our measurements of the sub-DLA at z ~ 2.2 to- 
wards HE 0001—2340 unveil a number of interesting abun- 



dance properties in this absorption system. With [0/H] = 
-1.81 ± 0.07, [N/H]< -3.3, and [N/0]< -1.5 the sub- 
DLA towards HE 0001— 2340 exhibits a classic massive 
star abundance pattern with little sign of primary ni- 
trogen (e.g., Umeda & Nomoto 2004). While the alpha- 
element abundances [O/H] and [Si/H] are comparable to 
those measured for other sub-DLAs and DLAs at similar 
redshifts (e.g., Dessauges-Zavadsky et al. 2002; Pettini et 
al. 2002), the [N/a] ([N/O]) ratio is among the lowest ever 
measured in DLAs and other astrophysical sites. [N/a] is 
a particularly important parameter to trace the chemical 
evolution in the early Universe. Nitrogen is believed to be 
produced mainly in intermediate-mass stars. The nitro- 
gen deficiency observed in the intergalactic medium and 
in blue compact dwarf (BCD) galaxies (e.g., Prochaska et 
al. 2002; Izotov & Thuan 1999) probably is a result from 
the early chemical enrichment of Type II supernovae from 
massive stars that produce only very little primary nitro- 
gen. In low-metallicity environments with [a/H] < —1.5, 
a first [N/a] plateau possibly is present near —1.5 dex, 
whereas for more metal-rich systems a second plateau oc- 
curs near —0.9 dex (Centurion et al. 2003). The region 
around [N/a] = —0.9 is much more populated than the 
region near —1.5 dex, showing that most of observed en- 
vironments must have experienced a second epoch of ni- 
trogen and a-element enrichment after the initial enrich- 
ment to the [N/a] = — 1.5 level. In contrast to the —0.9 
dex plateau, that near —1.5 dex exhibits only small scat- 
ter (~ 0.05 dex) in (N/O). To our knowledge, no systems 
with [N/O] ratios < —1.70 have been found so far. This 
may suggest that the value of [N/0]« —1.5 represents the 
ground floor for (N/O) in the Universe, set by the coeval 
injection of primary nitrogen and oxygen from the first 
generation of (massive) stars. The results presented here 
support this scenario. However, the number of systems 
with very low (N/O) values is still very limited and thus 
the existence of the —1.5 dex plateau still has to be con- 
firmed by additional measurements (Centurion et al. 2003; 
Israelian et al. 2004). 

Another interesting result of our study is the abnormal 
abundance pattern seen in components O and P. The ion 
column densities in these components derived in Sect. 3.4.2 
and the CLOUDY model presented in Sect. 4.2 imply sig- 
nificant overabundances of aluminum (0.5 dex), silicon 
(0.6 dex), and phosphorus (1.5 dex) compared to carbon. 
Such an abundance pattern is highly unusual for DLAs. 
The findings may indicate that the line of sight towards 
HE 0001 — 2340 passes through a confined gaseous region 
in the DLA host that is locally polluted with heavy ele- 
ments by nearby supernova explosions. Since the ejecta of 
a single supernova do not contain enough material to sig- 
nificantly enhance the heavy metal abundance in larger 
volumes of interstellar gas, it seems plausible that the 
chemically enriched gas components O and P is related to 
massive (coeval) star-formation activity, such as in mas- 
sive star clusters. Recent models for the so-called "hyper- 
novae" (supernovae that have explosion energies > 10 52 
erg) predict significantly enhanced [Si/C] abundances in 
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the supervova ejecta as a result of carbon and oxygen 
burning in low-density regions that lead to an enhanced 
freeze-out of a elements (e.g., Nakamura et al. 2001). Yet, 
these models cannot explain overabundances of phospho- 
rus and aluminum and thus the exact explosion mecha- 
nism that may lead to the observed abundance pattern in 
components O and P remains unknown. 

A striking feature of the z = 2.187 absorber towards 
HE 0001—2340 is the large number of individual velocity 
components that span a very large velocity range of al- 
most <~ 400 kms -1 . This velocity span seems very large 
for a single (proto)galactic structure. One possible expla- 
nation for this is that the absorption pattern reflects a 
merger event of two galaxies or protogalactic structures. 
The individual velocity components that arc quite narrow 
but that do not contain large columns of gas then may 
arise in tidal gaseous features produced by the merger 
event, spread out over a large spatial and radial- velocity 
range. An interesting aspect of a merger scenario lies in 
the idea that the merger actually triggers the formation 
of the stellar clusters (e.g., a nuclear cluster) that give 
rise to the abundance anomalies in components O and 
P. It is well known that merger/accretion events are of 
prime importance for starburst activity and the forma- 
tion of massive stellar clusters in and around galaxies 
(e.g., Knierman et al. 2003; Schweizer & Seitzer 1998). If 
the z = 2.187 absorber towards HE 0001 -2340 currently 
undergoes a merger event that leads to the formation of 
massive stellar clusters, these clusters should locally en- 
rich their interstellar environment with heavy elements by 
supernova explosions. Such circum-cluster environments, 
when they are young, should not have had enough time 
to mix with the more extended interstellar environment 
of the host galaxy. This possibly creates distinct regions 
with enhanced heavy-element abundances similar to what 
is observed in components O and P. The sub-DLA to- 
wards HE 0001— 2340 thus may represent a young galaxy 
at z w 2 that currently forms due to merging. 

6. Summary 

We have performed a detailed analysis of the sub-DLA at 
z = 2.187 towards HE 0001-2340 and have derived the 
following results: 

(1) We identify 16 absorption subcomponents in 
neutral and weakly ionized species (Hi, Di, Cn, Ni, 
Nil, Oi, Mgn, Alii, Aim, Sin, Pn, and Fen), and 
25 subcomponents in Civ and Siiv. These components 
span a velocity range of ~ 400 kms -1 in the z — 2.187 
restframe. The system has a total neutral hydrogen 
column density of log iV(Hi) « 19.7. 

(2) The overall abundance of the absorber is 
[M/H]= -1.81 ± 0.07 dex (1/65 solar), as derived 
from unsaturated O I absorption in the strongest neutral 
gas component A at v Ic i = +13 kms -1 (relative to 
z = 2.187). The lack of depiction of Sin and the absence 



of molecular hydrogen imply that the absorber is devoid 
of dust. 

(3) The gas in component A shows a classic mas- 
sive star abundance pattern with little sign of primary 
nitrogen. Taking possible ionization effects into account, 
we estimate [N/H]< -3.3 and [N/0]< -1.5. These values 
place the z — 2.187 absorber among systems with the 
lowest measured [N/a] ratios in astrophysical sites so far. 
This result is in line with the idea that the first stars in 
the Universe have enriched the intergalactic medium to a 
level of [N/a] w -1.5. 

(4) Components O and P near —350 kms -1 exhibit 
significant overabundances of phosphorus, silicon, 
and aluminum ([P/C]w +1.5, [Si/C]« +0.6, and 
[Al/C]« +0.5) after applying an ionization correction. 
Possibly, the gas belonging to these components is locally 
enriched by the supernova ejecta from one or more 
massive stellar clusters. Absorption by Di is possibly 
detected in components O and P in the blue wings of 
the Hi Lya, Ly/3 and Ly7 lines. However, our analysis 
of these lines leads to an extraordinarily high D/H ratio 
of > 1.22 x 10~ 4 , implying that the Di absorption is 
blended with an H I interloper. 

(5) We suggest that the large number of absorption 
components spread over several hundred kms -1 indicates 
the merging of two galaxies or protogalactic structures. 
This merging eventually triggers the formation of star 
clusters that are responsible for the abundance anomalies 
observed in components O and P. If so, the sub-DLA 
towards HE 0001— 2340 may represent a young galaxy at 
z « 2 that currently forms out of a merger event. 
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Appendix A: Comments on Deuterium 

The H I profile fit described in Sect. 3.3 (Fig. 5) suggests 
the presence of additional absorption in the blue wing of 
the Lya, Ly/3, and Ly7 profiles - exactly at the position 
where D I absorption from the two bluer-most components 
O and P is expected. A Voigt-profile fit to the possible D I 
Ly/3 and Ly7 absorption (shown in Fig. A.l) provides log 
A^(D l) = 13.94±0.04 and log N(D i) P = 13.54±0.05 to- 
gether with 6(D l) = 10.8 ± 1.7 and b(D i) P = 11.6 ± 1.9. 
The values for log iV(Di) do not change significantly if 
one assumes lower b values. With iV(Hi)op < 18 (see 
Sect. 3.3), we derive a limit for the deuterium-to-hydrogen 
ratio of (D/H) p > 1.22 x 1(T 4 . This is - 4 times 
higher than the canonical value for chemically young gas 
(3.0 x 10" 5 ; e.g., Buries et al. 2001). The simplest and 
most likely explanation for this unusually high (D/H) ra- 
tio is to assume that the deuterium lines in components O 
and P are blended with an 'hydrogen interloper', a weak 
Hi absorber that would be located near —440 kins" 1 (in 
the z = 2.187 restframe) and that coincidentally falls to- 
gether with the D I absorption in components O and P. 
However, as recently reviewed by Jedamzik (2002), there 
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Fig. A.l. D I Ly/3, and Ly7 absorption in the z = 2.187 
absorber. The solid line represents the best fit to the data 
with log AT(Di) = 13.94 ± 0.04 for component O and 
13.54 ± 0.05 for component P assuming b = 10.8 ± 1.7 
kms -1 and 11.6 ± 1.9 kms -1 , respectively. 

are a number of possible production mechanisms at high 
z that may locally enhance the (D /H) ratio in certain en- 
vironments, e.g., in regions predominantly enriched by su- 
pcrmassive stars. It is an intriguing coincidence that the 
apparently enhanced (D/H) ratio in the z = 2.187 ab- 
sorber toward HE 0001— 2340 is observed in gas that ex- 
hibits abundance anomalies also in other elements (P, Si, 
and Al; see Sect. 4.2). We therefore have to acknowledge 
that the existence of an intrinsically enhanced (D/H) ratio 
in components O and P cannot be entirely excluded, al- 
though the presence of an hydrogen interloper certainly re- 
mains the most likely explanation for the apparent (D /H) 
anomaly. 



